Effect of external electric field on the charge density waves in one 

dimensional Hubbard superlattices 
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We have studied the ground state of the one dimensional Hubbard superlattice structures with 
different unit cell sizes in the presence of electric field. Self consistent Hartree-Fock approximation 
calculation is done in the weak to intermediate interaction regime. Studying the charge gap at the 
Fermi level and the charge density structure factor, we get an idea how the charge modulation on 
the superlattice is governed by the competition between the electronic correlation and the external 
electric field. 

PACS numbers: 73.21. Cd, 71.30.+h, 71.45.Lr, 74.25. Jb 



I 

o 
o 



> 
in 



oo 

o 



I. INTRODUCTION 

The study of low dimensional metallic multilayered 
structures jlj is interesting because of their unique 
characteristics. The oscillation of exchange coupling 
between magnetic layers [2[ and the appearance of giant 
magnetoresistance [3| are among the exciting features 
of the multilayers. To investigate the properties of 
the metallic multilayers many theoretical works have 
been done taking simple superlattice structures as the 
models [3, 0, @, 13, [1] ■ This kind of models consist of 
periodic arrangement of Nu sites with repulsive on-site 
Coulomb interaction U (> 0), followed by A^o sites 
with no on-site interaction {U = 0). Some of these 
works investigated the ordering 0, 0| of the ground 
state, while some other explored the possibility of the 
metal- insulator transition in these systems 0j 9 ■ There 
are possibilities of formation of novel ground states such 
as charge ordered or spin ordered ones depending on 
the distribution of the interaction parameter U in such 
superlattices. 

On the other hand, the effect of electric field on the 
strongly correlated low dimensional electronic systems 
has attracted much interest in recent years because 
of their practical applications in tuning dielectric and 
piezoelectric properties Q. Many experiments are done 
on these low dimensional systems in the presence of 
electric field. It was found that spin ordered or charge 
ordered phases of a Mott insulator collapse in an electric 
field [H [U, m, m, [11 . Also some theoretical works are 
done on such systems in which a uniform electric field is 
implemented in the form of a ramp potential. Applying 
such an electric field in the homogeneous Hubbard model 
it was found that the field can induce oscillations in the 
charge gap of these systems [l5|, ll6|. However, it is not 
yet known how the superlattice systems behave in the 
presence of such an electric field. 

In this work, we investigate the electronic properties of 



simple superlattice structures in the presence of electric 
field. We consider weak to intermediate interaction 
regime and work under the Hartree-Fock Approximation 
(HFA). 



II. THE MODEL AND THE HARTREE-FOCK 
APPROXIMATION 

Our model is a one dimensional iV-site Hubbard chain 
with open boundaries. The model Hamiltonian is. 



H = ^ arii +^y^(4.o-Ci+i,, 



(1) 

where c\ ^{ci^a) is the creation(annihilation) operator 
for an electron with spin a (f or j) at the i-th site. 

fii,c7 = c|g.Ci,o-, and rii — X^o- "i.o' number opera- 

tor at the j-th site; t is the hopping integral between the 
nearest neighbor sites. Ui denotes the on-site Coulomb 
repulsion at the i-th site; in a superlattice UiS follow a 
repeated pattern depending on the size of a unit cell of 
the superlattice. is the site energy of the i-th site. In 
the absence of electric field all e^'s are set to zero. The 
external electric field is applied on the system in the form 
of a ramp potential J^] . In presence of this field the site 
energies become 
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where V is the applied voltage. This form of the site 
potential is used to ensure that the external bias varies 
from —V/2 to V/2 across the superlattice. We will work 
in the weak to intermediate coupling regime where U^t. 
It was observed in a previous work that in this regime 
the mean field approximation is quite reliable for this 
class of systems [3]. 



2 



We decouple the Hamiltonian using the unrestricted 
Hartree-Fock Approximation (HFA), 

(2) 

where (• • •) denotes the expectation value calculated 
with respect to the ground state. Now the Hamiltonian 
can be divided into two parts for the two types of spins, 
i.e. H = H-f + Hi. In an unrestricted Hartree-Fock 
Approximation, one determines the distribution of the 
rii^o-'s by diagonalizing and Hi in a self-consistent 
manner. The ground state is constructed by filling up 
the energy levels from both the up and the down bands 
upto the Fermi level. 




FIG. 1: Two different types of superlattice structures studied 
in the present work. The arrow shows the direction of the 
electric field. 

In this paper, we have presented the results for two 
different types of superlattice structures as shown in 
Fig. 1. For the first one the size of the unit cell is 
two and for the second one it is three. We have taken 
two types of sites with on-site correlation parameters 
U — 1 and U — respectively. We have studied the 
superlattices for various values of N. Since our aim is 
to compare the effects of the electric field on different 
superlattice structures of small size, we have presented 
here the results for the cases with N = 60 and 120 only. 
Comparing the results for these two system sizes one can 
also make an idea what happens in the infinite limit. 

To study the effect of the electric field on the metal- 
lic/insulating behavior of the ground state, we performed 
a systematic study of the charge gap (A) at the Fermi 
level of the system containing n electrons, 

A — En+l + En-1 ~ 2En, (3) 

where En is the ground state energy of an n electron 
system. We have also studied the charge density wave 
(CDW) structure factor 

hi 

where p is the average particle density on the super- 
lattice, Tp denotes the position of the p-th site and q 
is the wave vector. These two quantities enable us to 



capture the competition between the electric field and 
the correlation parameter in determining the charge 
modulation along the superlattice. 



III. RESULTS OF HFA CALCULATIONS 

It is well known that in the absence of electric field, a 
homogeneous Hubbard chain is an anti-ferromagnetically 
ordered system with a finite charge gap at half-filling 
[13, [III . Figure 2 shows the variation of the charge gap 
with the electric field for a half-filled homogeneous Hub- 
bard chain. It is clear from the diagram that the charge 
gap goes through a number of maxima and minima with 
increasing electric field. As N increases the minima shift 
towards lower values of V. The nature of variation of 
the charge gap is in good qualitative agreement with 
the previous Density Matrix Renormalization Group 
(DMRG) results [H. 




FIG. 2: Variation of the charge gap (A) at the Fermi level 
for the homogeneous Hubbard chain at half-filling with the 
applied voltage (V^). The solid line corresponds to A' = 60, 
while the dotted line shows the case of A = 120. 

Now we discuss the results for the —0 — U— super- 
lattice structure. At half-filling the system is a CDW 
insulator with g = tt in the absence of electric field [1, [3| 
and the charge gap is finite (V = case in Fig. 3(a)). 
As we turn on and gradually increase the electric field, 
the charge gap passes through a number of maxima and 
minima. For = 60 and 120thc variation of the charge 
gap are shown in Fig. 3(a). The type of oscillation 
of the charge gap with electric field is quite similar 
to that observed in the homogeneous system (Fig. 2). 
These oscillations are observable only in the finite sized 
systems and are crucially controlled by the interplay 
of the Hubbard interaction and the spatial gradient of 
the external bias. In the absence of electric field, the 
CDW phase (with q — n) at half-filling has tendency 
to form "doublons" (jj) at the sites with U — 0, while 
the sites with U > tend to depopulate. When the 
electric field is increased the electrons are pushed back 
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(a) -0-U- Model , U = 1.0 
1/2 filling 
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(b) -0-0-U- Model , U = 1.0 
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FIG. 3: Variations of the charge gap (A) at the Fermi level 
of the superlattice systems with the apphed voltage {V); (a) 
for —0 — U— superlattice at | filling, (b) for —0 — — U— 
superlattice at | filling and (c) for —0 — — U— superlattice at 
i filling respectively. The sohd line corresponds to N = 60, 
while the dotted line shows the case of N = 120. 



near one end of the lattice. This leads to population of 
some sites with U > 0, leading to the breakdown of the 
q = IT CDW phase. Such a crossover from the CDW 
phase is marked by the first minimum in the charge 
gap. Subsequent increase in the electric field results 
in gradual accumulation of electrons in one half of the 
superlattice. Because of the competition between the 
Coulomb correlation energy and the electric field term, 
the aforesaid process of piling up of electrons takes place 
only after finite increments of electric field in a finite 
sized system. This results in the oscillation in the charge 
gap. 

Next we study the CDW structure factor of this 
model. Generally speaking, from the mean field point 
of view, all sites with the same value of in a super- 
lattice are of same status. So a peak depending on the 
periodicity of the structure of the superlattice appears 
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(a) -0-U- Model , U=1.0 
V = 0.0 
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(b) -0-U- Model, U=1.0 
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FIG. 4: Charge density structure factor C{q) for the half filled 

— — U— superlattice systems; (a) for V — 0.0 and (b) for 
V = 1.25 respectively. Tiie solid line corresponds to N = 60, 
while the dotted line represents the case of N = 120. 



in the charge density structure factor. Our method 
also detects other peaks in C(q')'s that depend on the 
specific values of the density (i.e. the position of the 
peak depends on the Fermi wave vector kp)- As the 
system size N increases, the peak due to the structural 
periodicity of the lattice becomes larger compared to 
the other peaks. Under periodic boundary condition or 
in the limit of TV ^ oo, only the peak due to the struc- 
tural modulation survives and the other peaks disappear. 

Referring to the specific case of —0 — U— superlattice 
at half-filhng, there is a sharp peak in C{q) a,t q ~ tt in 
the absence of electric field (Fig. 4(a)). Incidentally, in 
this case the peak due to the periodicity of the lattice 
structure and the 2fcf peak both occur at (/ ^ tt. As 
the electric field is increased, after a critical value of 
the same, a peak at ~ is seen (see Fig. 4(6)). This 
is due to accumulation of the charges near one side of 
the chain. At this value of the electric field, the peak 
at q ^ n starts to get diminished. For larger electric 
field, the peak at g ~ becomes larger at the cost of 
the peak at g ~ tt. In Fig. 5(a) the variations of the 
magnitudes of the peaks at g ~ and q ~ tt are shown. 
It clearly shows that the value of the electric field at 
which the q ^ peak becomes significant is as same as 
the value of the electric field where the first minimum 
of the charge gap oscillation occurs (see Fig. 3(a)). 
So it is clear that at this value of the electric field 
charge accumulation at one side of the chain begins to 
dominate and the charge ordering tends to get destroyed. 
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Then the distribution of the site potentials in a unit cell 
of the superlattice is in unison with the distribution of 
the correlation parameter. This phenomenon reinforces 
the aforesaid charge ordering for small values of the 
electric field. So an increase in the value of the charge 
gap is observed. For larger electric field, however the 
spatial gradient of the site potentials becomes so large, 
that a global shift of the charges towards left is preferred 
in the chain and the charge ordering is destroyed. Also 
for the —0 — U— model a trace of such feature is observed 
for low electric field. There is a change in the slope of 
the charge gap as a function of the applied voltage, for 
a low value of the voltage (see Fig. 3(a)). Since in this 
case the unit cell contains only two sites, the feature 
is only weakly showing up. On the contrary, we have 
noted that for —U — — 0— model and —U — 0— model, 
this type of enhancement of charge ordering for small 
values of the electric field is absent. In those cases, on 
application of the electric field the charge gaps initially 
fall and then subsequent oscillations are observed. 



FIG. 5: Variations of the height of the peaks in C{q) for the 

superlattice systems with the applied voltage {V); dashed line 
is for tiic peak at g ~ 0, while solid line is for the peak at g ~ vr 
and g ~ ^ in (a) —0 — U— superlattice at | filling and (b) 
—0 — — U— superlattice at | filling respectively. 

This point can be taken as a transition point, though 
on either sides of it the system remains insulating. 
By observing the increase in height of the peak at 
g' ~ 0, one can easily understand how the charge accu- 
mulation at one side of the chain grows with electric field. 

In Fig. 3(6) we present the variation of the charge 
gap with the electric field for the other system, the 
— — — U— superlattice at 2/3 filling. At zero electric 
field, the electrons try to accumulate at the sites with 
zero Hubbard interaction as it minimizes the energy. As 
a result, the system is charge ordered insulator with a 
finite charge gap. A sharp peak of C{q) at g ~ ^ is 
obtained. Here also the density dependent 2kp peak 
and the peak due to the structural periodicity of the 
lattice occur at the same q (~ ^). On application of the 
electric field, the charge gap initially increases and then 
falls gradually. Subsequent oscillations are observed as 
in the previous case of —0 — U— model. 

The initial increase in the charge gap for small values 
of the electric field can be understood in the following 
way. In a particular cell of this —0 — — U— superlattice 
at 2/3 filling, the two sites at the left (with U = 0) are 
preferred by the electrons in the absence of electric field; 
they tend to be doubly occupied keeping the other site 
empty. In this situation a charge density wave is formed 
in the chain. For nonzero electric field, there is a positive 
gradient of the site potentials towards right (see Fig. 1). 
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FIG. 6: Variations of the chaxge gap (A) at the Fermi level for 
the quarter filled superlattice systems with the applied voltage 
{V); (a) for -0 - U- superlattice and (b) for -0 - - ?7- 
superlattice respectively. The solid line corresponds to N = 
60, while the dotted line shows the case oi N = 120. 

Figure 5(6) shows the variation of the heights of the 
peaks of C(g) with external bias for the —0 — — U— 
model; type of variation is quite similar to the previous 
case of —0 — U— model. The peak at g ~ ^ starts to 
fall and a peak at g ~ appears at a critical value of the 
electric field. This critical value again matches with the 
value of the electric field at which the first minimum of 
the charge gap oscillation occurs in Fig. 3(6). So in this 
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FIG. 7: Charge density structure factor C{q) for the quar- 
ter filled suporlattico systems at zero electric field; (a) for 
—0 — U— superlattice and (b) for —0 — — U— superlattice 
respectively. The solid line corresponds to N = 60, while the 
dotted line shows the case oi N = 120. 



case also it makes a transition from a CDW phase to a 
phase where electrons pile up at one end. 

We have shown the variation of the charge gap for the 
— — — [/— superlattice at half-filHng in Fig. 3(c) for the 
sake of comparison with the other half-filled cases shown 
in Fig. 2 and Fig. 3(a). Here at ^ = the gap diminishes 
with larger values of N indicating a metallic behavior at 
the thermodynamic limit. An increase in V drives the 
system ultimately to an insulating one. 

We have also studied the quarter filled cases for our 
superlattice models. Figures 6(a) and 6(6) show the 
charge gap oscillations in —0 — U— and —0 — — U— 
models respectively. Both the systems are metallic in 
absence of the electric field; the charge gap A ~ -i-, 
N being the system size. With the increase of electric 
field the charge gap oscillatc^s between a number of 
maxima and minima, and then increases monotonically, 
indicating an insulating phase. 

In Fig. 7 we plot the C(g)'s for both types of the 
superlattices at quarter filling in absence of the electric 
field. For —0 — U— model (Fig. 7(a)), we see a dominant 
peak at q = TT, which is due to the structural periodicity 
of the lattice. Apart from this, a 2kp peak appears at 
g ~ J. Also few other wiggles due to various possible 
short length scale density modulations are observed. 
For —0 — — U— model (Fig. 7(6)) also, the dominant 
peak is due to the lattice structure and it occurs at 



FIG. 8: Variations of the height of peaks in C{q) for the 

quarter filled superlattice systems with the applied voltage 
{V); dashed line is for the peak at g ~ 0, small dashed line 
is for the peak a± q ~ ^ (2k f) and solid line is for the peak 
due to structure (a) for —0 — f/— superlattice and (b) for 
—0 — — U— superlattice respectively. 



q ~ Here also other wiggles due to different possible 
density modulations are obtained along with the 2kF 
peak {q ~ |). 

It can be seen explicitly, in the limit U = 0, that 

the position of the 2kp peak in finite chains is shifted 
from the expected position q = 2kp = mT/{N + 1) (n 
being the number of occupied single particle levels) by 
an amount ~ As the system size N increases, they 
approach the actual 2kF value. On the other hand, the 
peak heights fall as So in the large N limit (and also 
under periodic boundary condition) the 2kF peaks are 
not detectable. 

As the electric field is turned on and increased, a 
peak at g ~ appears for both the models at quarter 
filling. The peak due to the structural periodicity and 
the 2kF peak (also the wiggles in C{q)) fall off; rate 
of fall is much slower for the structural periodicity 
peak. These show the gradual accumulation of the 
charges near one side of the chain. In Fig. 8 variations 
of the heights of the peaks in C{qys with the external 
bias is shown. It clearly shows how charge accumu- 
lation at one side of the chain increases with electric field. 
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IV. CONCLUSION 

In this work, we have studied the one dimensional Hub- 
bard superlattices with different types of unit ceUs in the 
presence of electric field. Single orbital nearest neigh- 
bor tight-binding model has been used. We have main- 
tained the fixed system sizes {N — 60 and 120), for the 
purpose of comparison of these models. The charge gap 
and the CDW structure factor of the systems are stud- 
ied under Hartree-Fock approximation in the presence 
of electric field. To check the reliability of the Hartree- 
Fock approximation results, we compared our results for 
the homogeneous Hubbard model with a previous DMRG 
calculation[15] and found reasonable qualitative agree- 
ment. Oscillations in charge gap obtained in the super- 
lattice systems are rather similar to those observed in the 
homogeneous Hubbard chains, showing signature of the 
finite sizes. Variations of the heights of different peaks 
of C{q) with the applied electric field give an idea about 
the distribution of electrons on the lattice. We found 
that the 2fci? peaks of C{q) arise for finite sized systems 
only. On application of the electric field a peak at g ~ 
appears at a critical value of the field and then increases 
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